clear

t = 0:0.01:30;

x0 = [0; 0.5];

x = lsode("mySys", x0, t);

size(x);

pos = x(:,1);
vel = x(:,2);

figure(1);
plot(t, pos);
xlabel("Time, t (sec)");
ylabel("Position, x (m)");
title("Position vs. Time");

figure(2);
plot(t,vel);
xlabel("Time, t (sec)");
ylabel("Velocity, v (m/s)");
title("Velocity vs. Time");
